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Abstract 


We study the kinetic regime of the Bose-condensation of scalar particles with 
weak self-interaction. The Boltzmann equation is solved numerically. We 
consider two kinetic stages. At the first stage the condensate is still absent 
but there is a nonzero inflow of particles towards p = 0 and the distribution 
function at p = 0 grows from finite values to infinity in a finite time. We 
observe a profound similarity between Bose-condensation and Kolmogorov 
turbulence. At the second stage there are two components, the condensate 
and particles, reaching their equilibrium values. We show that the evolution 
in both stages proceeds in a self-similar way and find the time needed for 
condensation. We do not consider a phase transition from the first stage 
to the second. Condensation of self-interacting bosons is compared to the 
condensation driven by interaction with a cold gas of fermions; the latter 
turns out to be self-similar too. Exploiting the self-similarity we obtain a 
number of analytical results in all cases. 
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I. INTRODUCTION 


It is a fundamental result of quantum statistics that above a certain critical density all 
added bosons must enter the ground state: a Bose-Einstein condensate (BEC) forms. Gas 
of weakly interacting bosons allows a well-dehned theoretical microscopic treatment of this 
process. However, Bose condensation has not yet been observed in a real physical system 
which could be described as an ideal gas. Superfluid "^He, for example, is a system with very 
strong interaction and is considered as a quantum liquid. 

One can reach Bose condensation by gradually decreasing the temperature in a sequence 
of equilibrium states. An appropriate description of this regime is given by the well-known 
theory of second order phase transitions []^. On the other hand, the conditions for the 
formation of a Bose condensate can appear when the system is far from equilibrium. This will 
be generally true for the condensation of an ideal gas because of the weakness of interaction. 
The kinetics of Bose-condensation in this case is an important and interesting problem. 

Recently this subject has attracted particular interest in connection with the exciting 
prospects for experimental observation of a Bose condensation in very cold atomic samples, 
e.g. in a gas of spin-polarized atomic hydrogen or in alkali-metal vapors as well 

as in the gas of excitons (bound states of electrons and holes in semiconductors) p. For 
instance, Doyle et. al. 0 report to have attained n ^ 10^^ cm“^ and T 100 /iK, which is 
only a factor of 3 above the temperature of Bose condensation at this density. 

Since the life-time of a sample in all experiments is hnite and short, the question of a 
time-scale for the formation of Bose-condensate is a central one. 

Another interesting application of Bose kinetics is rather far from laboratory experiments 
and is related to the problem of Bose star formation from the dark matter in the 
Universe. The life-time of the “sample” is not limited in this case, but in the only model 
for Bose-star formation which involves a realistic dark matter particle candidate, the axion 
IP], the self-interaction of particles is so weak that the relaxation time is comparable to the 
age of the Universe. The possibility of Bose-star formation becomes dependent upon the 
exact value of this time-scale. 


In experiments with atomic hydrogen or with alkali-metal vapors the Bose-gas is effec¬ 
tively isolated from its surroundings and the relaxation to thermal equilibrium and con¬ 
densation is expected to occur due to particles’ self-coupling only. The same is true for 
the case of axion miniclusters. On the contrary, interaction with the thermal heat bath is 
important in the dynamics of excitons. We consider both situations in the present paper. 
However, we only consider the simplest self-interaction of structureless scalar particles 
(this approximation is exact e.g. in the axion case). The heat bath is modelled by a gas of 
heavy fermions. 


The process of Bose condensation can be divided into three stages. The hrst and the 
third ones are kinetic stages which occur before and after actual nucleation of the condensate. 
These two stages can be treated with the Boltzmann equation which describes restructuring 
of the distribution function at the hrst stage and growth of the condensate at the third. A 
numerical study of these two stages is the purpose of the present paper (we have reported 
the main results in Ref. 
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In the framework of the kinetic equation there is no condensate at all times if there was 


2 





no condensate initially. Therefore we have to add a small seed condensate “by hand” in 
this framework when switching from the hrst to the third stages. We switch between stages 
when the distribution function becomes inhnite at zero momentum. We do not consider the 
second stage, i.e. how the condensate actually emerges, which is, in fact, a phase transition. 

ni. Since the duration of 
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A description of the second stage has been elaborated in Ref. 
this stage is short and the magnitude of the condensate which appears is small 
approach seems to be reasonable. 

The question of the time evolution of a weakly interacting Bose gas during the hrst 
kinetic stage was addressed in a number of papers. In earlier treatments an ideal Bose gas 
was coupled to a thermal bath with inhnite heat capacity: to a phonon bath in Ref. |jT^ and 
to a fermi heat bath in Ref. |TB]. A small energy exchange was assumed and after several 
other approximations a Fokker-Planck type of equation was obtained. Levich and Yahkot 
calculated analytically that the time for condensation (more precisely, the time needed 


for the distribution function to reach inhnity at zero momentum) is inhnite in this situation. 
By including boson-boson interactions they later [|T^] found a solution which describes the 
explosive emergence of a condensate, but they concluded that this ehect could have been an 
artifact of their approximations. 

Snoke and Wolfe in Ref. undertook a direct numerical integration of the Boltzmann 
kinetic equation. Although this calculation demonstrated a restructuring of the distribution 
of particles, the appearance of a Bose condensate (more precisely, an inhnite value of the 
distribution function at zero momentum) was not detected. As compared with Ref. JT^, we 


carry out the numerical integration in a much wider dynamical range of relative energies and 
densities. The major and important diherence of our approach is that we directly analyse 
the behavior of the distribution function, while Snoke and Wolfe considered an integrated 
quantity n{e), the number of particles with the energy less than a given value e. 

Analytical study of Bose condensation was performed recently by Kagan, Svistunov and 
Shlyapnikov in the paper |]R| where three diherent regimes of evolution were identihed and 
considered. Our numerical results conhrm qualitative predictions of that work. In particular, 
it was argued that in the kinetic region of the non-linear regime the distribution function 
has the form f{e) ~ for small e. We indeed observe the tendency to this law in our 
numerical integration, but, while being close to it, the system never reaches this critical 
exponent. Instead of 7/6 ~ 1.17 we obtain 1.24. Moreover, this distribution is destroyed 
when the condensate emerges, the fact that was not stressed in Ref. 


161. As concerns the 


condensation time, only order of magnitude dimensional estimates were obtained in Refs. 

00 . 

We hnd that during the hrst kinetic stage the distribution function became inhnite at 
zero momentum in a hnite time in the case of self-interacting bosons. This stage can be 
divided into two parts. During the hrst part of the hrst stage the distribution function is 
restructured so that a self-similar solution is formed; the distribution function does not grow 
much during this time interval which takes about half of the entire evolution time. During 
the second part, the distribution function rapidly (as a pole power law oc — reaches 

inhnity at p = 0 in a self-similar way. 

Despite the fact that the Bose-condensation is inherently based on Bose-statistics, the 
kinetics of this stage is a classical process and not a quantum one. To be precise, it corre¬ 
sponds to the kinetics of classical waves (but not classical particles). Actually, the equations 
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we encounter here and the resulting behavior of the system are well-known in the theory of 
turbulence ||T7|,^,^ . This is no wonder. Indeed, a Bose condensation occurs when the occu¬ 


pation numbers are large, /(p) ^ 1. In that case we can replace quantum creation operators 
by amplitudes of classical waves. The corresponding kinetic equation can be obtained in 
two ways: either neglecting 1 in the product of all factors [1 -|- /(pi)] in the quantum kinetic 
equation (by which we neglect the spontaneous scattering with respect to the induced one). 


or directly in the random phase approximation for the ensemble of classical waves |21 


The second approach is customary in the theory of turbulence, but the result is the same 
in both approaches.^] While it may be more appropriate in the second approach to refer to 
f{e) as a density of waves in phase space, we will continue to call it density of particles in 
phase space. 

The power law behavior of the distribution function we observe, with an exponent close to 
7/6, is a special case of the Kolmogorov turbulence il8hp0|. Indeed, by dehnition, turbulence 
is a stationary state characterized by a non-zero flow (in phase space) of some conserved 
quantity. In the case of Bose-condensation we have a flux of particles towards the condensate, 
which tends to be constant across the energy space, and this corresponds to f{e) oc 
An explosive behavior of the distribution function at p = 0 which we observe was also found 


in other instances of turbulence 22,23 . 


The actual nucleation of the condensate cannot be described in the framework of the 
kinetic equation. But we have to conclude that after the distribution function has become 
inhnite at zero momentum, the condensate forms. Kinetic equations which describe the 
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We re-derive 


third stage in the presence of the condensate were derived in Refs, 
those equations in a simpler way which is valid when the difference between particles and 
quasiparticles is unimportant. We consider this weak coupling regime only. We found that 
during the third stage the evolution is also self-similar until the restructuring in the phase- 
space reaches the tail of the distribution. 

We compare this condensation process, which is due to self-interaction of particles, to 
the condensation driven by the interaction with a cold bath, for which we also solved the 
Boltzmann equation numerically. We hud that the distribution function grows with time only 
as oc in this case. The energy dependence of the distribution function is also different now. 
Instead of a singular power law prohle, f{e) ~ which corresponded to the constant 

flux of particles towards the condensate, now all excess particles form a well dehned packet in 
the momentum space. The total number of particles in this packet remains constant during 
the evolution, while its width decreases, gradually approaching a 5-function [^. This type 
of behavior is qualitatively different from the one we have in the case of self-interacting 
bosons. However, it is also self-similar and one of the predictions of the self-similarity in this 
case is that f{e = 0,t)(xt^. In general, using the properties of the self-similarity enables 
us to obtain a number of analytical results, e.g. the relation of the critical exponent (3 in 
f{e = 0, f) oc {tc — t)~^ to the critical exponent a in f{e) oc e~°‘. 

The fact that the system interacting with a cold bath needs inhnite time (in the kinetic 


^This is particularly important for the axion case, because axions in miniclusters [0 are from the very 
beginning in a state of gravitationally bounded classical waves and one may wonder (or be suspicious of) 
how Bose-condensation could occur in such a system. 
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framework) to reach infinity in the distribution function at p = 0 does not necessarily 
mean that the condensation time is inhnite for this system. Indeed, the nucleation of the 
condensate is a phase transition which does not have to occur at the moment when the 
distribution reaches inhnity, and can occur earlier. 


Although our prime interest and motivation for this work were connected with the physics 
of Bose stars, the problem of condensate formation in a gravitational held has never been 
considered, and we do not consider it here. Instead, we are solving kinetic equations in hat 
space-time and only the range of parameters of the initial distribution can rehect the virial 
equilibrium of a self-gravitating system. 


When this work was completed B. Svistunov brought to our attention his paper |^. He 
also analyzed two non-linear kinetic stages exploiting the assumption of self-similarity. His 
approach is close to ours and our results qualitatively conhrm the results and predictions of 
Ref. However, he assumed a priory in his analysis that the critical exponent is a = 7/6. 
He concluded that this assumption had passed his numerical self-consistency check, but 
we believe that this is due to the fact that he used energy grid with the minimum value 
£min = 10“^, while in our case it was Emm = 10“®. 


The plan of this work is as follows. In Sec. y we reduce the Boltzmann kinetic equation 
to the form suitable for numerical integration, both in the absence and presence of the 
condensate. In Sec. El we present the results of numerical integration which describe the 
evolution of one-particle distribution function in the case of self-interacting bosons, while 
Sec. 1^ is devoted to the evolution of a Bose system interacting with fermion bath of inhnite 
heat capacity. In Sec. |V| we derive analytical consequences of self-similar behavior of the 
distribution function. In the hnal section 0 we discuss some possible applications and 
present conclusions. Appendix y describes an analytical reduction of the collision integral. 
In Appendix y the static (“turbulent”) power law solutions to the kinetic equation are 
described. Appendix P is a technical supplement to Sec. |V. 


II. THE KINETIC EQUATION 


In a state of thermal equilibrium any system of bosons has to have the Bose-Einstein 
distribution of particles over the energies. In order to reach this state particles need to 
interact either with the thermal bath of other particles or with each other. We shall consider 
both cases: the system of scalar bosons with 4-particle self-interaction and an ideal Bose 
gas coupled to cold gas of fermions. The time development of a state which contains large 
number of particles can be adequately described by the kinetic equation which governs the 
evolution of the one-particle distribution function, /(p). The only process which contributes 
to the Boltzmann kinetic equation is illustrated in Fig. for each case we consider, and the 
equation takes the form: 


#(Pi) 

dt 


2poi 




k 


k=2 


d^Pk 

(27r)32pofc ’ 


( 2 . 1 ) 


where {pok,Pk) are components of the four-vector p^k- For the conceptual convenience we 
wrote this equation in the most general Lorentz-invariant form, but in what follows we shall 
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consider only the non-relativistic limit of it. For the qnartic self-interaction, Fig. |I|a, we 
have 

^if) = [1 + /l][l + / 2 ]/i /2 “ [1 + /l][l + /2 ]/i/ 25 (2.2a) 

where fk = f{Pk), fk = /(PfcO assnme the convention 1' = 3, 2' = 4), while for the 
interaction with the fermion bath, Fig lb : 

F{f) = [1 + /i][l - V2]fW2 - [1 + /{][1 - ^2]/i^2, (2-2b) 


where % = ri{pk) is the one particle distribntion fnnction of fermions which we shall not 
evolve here and consider it to be fixed to the eqnilibrinm. 

The field-theoretical Lagrangian which we bear in mind for the case of self-interacting 
bosons is 


C = 



(2.3) 


and the corresponding matrix element is given by 

= A". 


(2.4) 


The matrix element of the process in Fig. lb in the non-relativistic limit is also momentnm 
independent, = g‘^. 

The Bose-Einstein distribntion fnnction is 


/(P) 


1 

exp [(£ - fi)/T] - 1 


(27r)V(5(p) , 


(2.5) 


where e is the particle energy, /i is the chemical potential, T is the temperatnre of the hnal 
state and ric is the density of particles in the condensate. The fact that distribntion fnnction 
( p.5|) is a solntion to the eqnation df{pi)/dt = /coii(Pi) = 0 is a simple conseqnence of 
the energy-momentnm conservation in each particle interaction vertex and does not depend 
npon the specihc form of the matrix element Mfi. 


One expects that an arbitrary initial distribntion fnnction /(p) will evolve via eqnation 
( p.l|) to its static solntion ( p.5|) just because the static solution can not evolve any further. [| 
This implies that simple and straightforward numerical iteration procedure has to be robust. 
Given the distribution function at time t we calculate the collision integral and then we 
advance the distribution function in time as /(p, t + At) = /(p, t) + /coii(p, t)At. 


The evolution of the initial distribution function which allows for the condensate for¬ 
mation consists of two stages. In the first stage condensate is absent and the distribution 
function is hnite. In the second stage ric ^ 0 and the distribution function has h-function 
singularity. Consequently, we have to derive the kinetic equation in the form suitable for 
the numerical integration in each of those cases separately. 


^The situation is not that simple, however, since the kinetic equation at / ^ 1 has many static solutions, 
and our system first tend to evolve to the turbulent distribution, see Appendix and would stay there 
(which it does in our numerical experiments) if the phase transition with a nucleation of the condensate 
would not occur. 
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A. Zero initial condensate 


In what follows we shall consider an isotropic initial distribution f = f{e) and the non- 
relativistic limit £ = p^/2m. The kinetic equation for the case without condensate, Uc = 0, 
can be rewritten in the form (see Appendix ^): 


dfjsi) 

dt 64:71^771 


F{f)-de[de', ^ Ip , 

Pi 


( 2 . 6 ) 


where 


D = min[pi,p2,pi,P2] • (2-7) 

The integration should be done over the white area of Fig. ^ which is divided into four 
regions. In each of those regions we have 


D = P 2 in I region 

D = Pi in II region 

D = P 2 in III region 

D = p'l in IV region 


In the argument of F{f) according to the conservation of energy 82 = + £2 ~ 

It is possible to map regions II-IV in Fig. ^ to the region I , and then to rewrite 
the collision integral in Eq. (|2.6|) as an integral over the region I only. This representation 
is essential in Ending the stationary solutions to Eq. m 0 and the corresponding 
Kolmogorov’s spectra, see Appendix |BT. 


B. Non-zero condensate 


This case we shall consider for bosons with self-interaction only. After the moment of 


condensate formation the kinetic equation (|2.6|) is inappropriate for numerical integration, 
and the finite number of particles in the condensate corresponds to the infinite value of the 
distribution function at zero energy. In order to describe the system of particles interacting 
with the condensate we divide the distribution function into two pieces: 


/ = /(£,«) + (2ir)V(i)£(p) 


( 2 , 8 ) 


where the first term corresponds to the ”gas” of particles and the second one describes the 


condensate. Substituting this function into the original kinetic equation (p.l|) we obtain 


ndt) = 


A^ne(t) r 

64 : 7 T^m Jo 


de[de2[f[f2 — /2(l + /i + /2)] 


(2.9a) 


/(^i) = Ip + 


ne(t)A^ 
327 im‘^pi \Jo 


r£l roo 

/ [/1/2 “ /l(l + /l + /2)]'^^2 + 2 / [/2(1 + /i + / 2 ) — /l/2]d£2 

JO 2ei 


(2.9b) 


7 













In Eq. ( p.9b|) e[ = ei — e '2 and £2 = ^'2 ~ while in Eq. (p.9a|) 82 = 8 \ + e^- 

In general, after the condensate formation (and at large particle densities even before) the 
description in terms of qnasiparticles rather then particles is more appropriate. For example, 
the kinetic eqnation ( ^.91 ) does not include the processes where one of the incoming and one 
of the outgoing particles have zero momentum, P 2 = P 2 = 0. This process does not contribute 
to the collision integral directly; he., it does not change the distribution of particles over 
energies, but it does change the effective particle mass. However, in many cases those effects 
are insignihcant and we still can work in terms of particles. Quantitatively, the original 
particle picture is correct if the potential energy contribution from the interaction with the 
condensate to the effective particle mass squared, = A0^/6 = An/m, is smaller than 

the particle momentum squared, he. An/m In fact, the kinetic equations of Refs, 

which contains all quasiparticle effects do coincide in this limit with our Eq. (0) derived 
in a simple way. Since n ~ m^{Av)^fo, we obtain Xfo{Av) -C 1 as a condition of validity 
of our equations, where An is the characteristic velocity dispersion. In the case of axion 
miniclusters, for example, we have |§] A/o ~ 10“^; An ~ 10“®, and the description in terms 
of particles is perfectly good. 

Note also that the kinetic equation itself becomes invalid at characteristic energies smaller 
than the inverse relaxation time, but this condition is always weaker than the previous one. 




III. SELF-INTERACTING BOSONS 


A. Initial distribution and the rescaling of time 


Let us hrst consider the system of self-interacting bosons. As an initial distribution we 
choose the function f{e), which has a maximum at p = 0. In general, such a distribution 
function can be characterized by means of three major parameters: 

(1) The overall amplitude /q. In what follows we dehne /o = f{e = 0). 

(2) The energy scale eo, where the distribution function becomes 2 times smaller, /(^o) = 
/0/2. 

(3) The effective width T of the region over which the distribution function varies rapidly. 
More specihcally, we choose the initial distribution function to be of the form: 

f{e) = arctan [exp (r(l — e/eq))] . (3.1) 

TT 


In what follows we shall measure the distribution function in units of /o, he. the initial 
distribution function which will appear throughout the rest of this section will have the 
normalization /(e = 0) = 1, and we shall measure the energy in units of Eq. 

We dehne the dimensionless time r as 


go/o(l + /o)A^ 
bdTT^m 


(3.2) 
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The parameter Eq after rescaling will not enter the kinetic equation explicitly, but will 
dehne the time scale of the evolution. The parameter /o dehnes the time scale as well, but 
it remains present in the equation through the dehnition of F if /q ~ 1: 


FU) 


f'Ji - hk + /o([/i + k]![!', - I/; + /a/1/2) 

1 + /0 


(3.3) 


Otherwise it also disappears from the rescaled equation: in the limit /q -C 1 one can neglect 
terms ~ /^, and in the limit /o S> 1 it is possible to drop terms ~ p in Eq.(|3.3|) 

In terms of r there remains the weak dependence of the relaxation time upon the initial 
shape parameter T. All data presented in this section will correspond to one and the same 
value of T = 5. 

After rescaling everything in Eq. (p.6|) is of order unity and, very roughly, it is expected 
that the time-scale of the relaxation towards the equilibrium corresponds to r ~ 1, see Ref. 

0 . 

While the distribution function Eq. (|2.5| ) has three characteristic parameters {T, p, nd, 
only two of them are non-zero simultaneously. Those two non-zero parameters are uniquely 
dehned by the values of two conserved quantities: densities of the energy and the number 
of particles: 


e = 


n = 



(3.4a) 

(3.4b) 


With the shape of the initial distribution function being given, the two scaling parameters 
/o and So in addition to hxing the time scale dehne also the parameters of hnal equilibrium 
This correspondence is illustrated in Fig. 

The solid curve divide the plane of Fig. in two regions and corresponds to the equi¬ 
librium distribution function with p = nc = 0. Below this curve the chemical potential /i 
remains zero, but ric > 0. /.e., if the initial number density of particles and the energy 
density correspond to this region, then the condensate forms in the hnal state. On the 
contrary, for the region above the solid curve, p < 0 and ric = 0 in Eq.( p.5|) . The dashed 
line on the Fig. ^ corresponds to the variation of the parameter /q (keeping the parameter F 
of the initial distribution (pT|) to be hxed, F = 5). With the increase of /o the phase point 
moves along the dashed line from left to right. The choice /o ^ /crit ~ 2.8 corresponds to a 
Bose-gas density less than critical, while /o > /cut corresponds to condensate formation in 
the hnal equilibrium state. 


B. Kinetics without condensate formation in the equilibrium state 

For completeness (and as a particular check of our numerical procedures) let us hrst 
consider the case /o ^ /crit- We integrated kinetic equation with two particular values of 
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/o, /o = 0.1 and /o = 1. In the case /o = 0.1 the p qnantnm terms in the fnnction F{f), 
Eq. (|3.3|) , are small corrections and the kinetic eqnation rednces to the Boltzmann eqnation 
for classical particles, while /o = 1 is a trne qnantnm regime. 


The resnlts of nnmerical integration of the kinetic eqnation are presented in Fig. 

The dashed cnrve corresponds to the initial distribntion fnnction (|3.1|) . The dotted cnrves 
represent nnmerical resnlts at a seqnence of moments r 1,5,10,20. The equilibrinm 
distribntion is shown by the solid cnrve. It is close to the classical Boltzmann distribntion 
in Fig. and corresponds to the Bose-Finstein distribntion ( p.5|) with T 0.63 £o cind 
p —0.078 £o in Fig- |b. We see that already at r ~ 10 the nnmerical resnlts are almost 
indistingnishable from the eqnilibrinm distribntions, and by the time r ~ 20 they coincide 
identically on the scale of Fig. 0. 


On Fig. 1^ the valne of the distribntion fnnction at p = 0 is shown as a fnnction of time 
for both cases /o = 0.1 and /o = 1. We see that the relaxation time is close to 10 in the 
hrst case, while ~ 20 in the second case. However, according to Fq. (|3.2|) the physical 
relaxation time in the case /o = 1 is an order of magnitnde shorter, which is the effect of 
the stimnlated scattering. 


C. Non-linear quantum kinetic regime or the classical turbulence 


Now we tnrn to the case we are interested in: the Bose-condensate formation. The initial 
distribntion have to correspond to a large valnes of the parameter /o, e. g. in onr case of 
Fq. (B /o > 2.8. We shall simplify the problem and consider /o S> 1. In this case we can 
disregard terms in the fnnction (|2.2a|) , which becomes: 

F(f) = i/i+ fiiAA - [/:+/a/iA. (3.5) 


Parameter /q disappears from the kinetic eqnation and the latter became scale invariant. 
This scale invariance gives the origin to a self-similar evolntion, as we shall see shortly. 


As we have noted in the Introdnction, very large occnpation nnmbers corresponds to 
the statistics of classical waves. The kinetic eqnation with F{f) given by Fq. ( |3.5|) can be 
obtained directly in a random phase approximation for the classical waves and is well known 
in the theory of tnrbulence 


17,21 


We integrated the kinetic eqnation in the energy interval 10“® < £ < 10. We had dehned 
the distribntion fnnction on the grid of 200 points eqnally spaced in the logarithm of energy 
and nsed the spline interpolation when calcnlating the distribntion fnnction at intermediate 
points. We have checked that the grid of 400 points prodnces essentially identical results. 
For each integration in collision integral we had implied Gauss algorithm. Particle and 
energy non-conservation was of order 10“^ for the entire time of integration. 


The results of the numerical integration of the kinetic equation (|2.6| ) are presented in 
Fig. I, where we plot the distribution function at different moments of time. We arranged 
the output each time that f{e = £min, t) had increased by one order of magnitude. The most 
striking feature of this plot is the self-similar character of the evolution. The distribution 
function has a ’’core” where f{e) ~ const and the radius of the core decreases with time 
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while the value of f{e) in the origin grows. Self-similar solutions exhibiting this kind of 
behavior can be parametrized as f{e,T) = A~^{T)fs{e/A{T)), where it is assumed (3 > 0 
and A{t) —>■ 0 with the increase of time. Outside the core the distribution function does not 
depend upon time and is the power law to very good accuracy. 

The value of the logarithmic derivative of f{e) is plotted in Fig. The dotted lines in 
Figs. 6 correspond to the limiting value a = 7/6 (see Appendix but this value is never 
reached prior to the moment of condensate formation. After that moment the character 
of the evolution completely changes (see the next Section). Rather, with the boundary 
condition df /de = 0 at e = 0 the power law on the tail is a = 1.24. 

We believe that the difference of the exponent in our self-similar distribution from the 
limiting critical exponent a = 7/6 is not a kind of the numerical artifact. In favour of that 
we have the following arguments: 

1. Doubling the number of grid points and/or integrator accuracy does not change the 
distribution function appreciably. We have tried also different integration algorithms and 
different boundary conditions at large energies. 


2. We had observed that after the particle flux wave reaches the bottom boundary of 
the integration interval emin = 10“® (which happens at r ~ 19 in our case), the distribution 
rapidly relaxes to the strict power law throughout the whole integration interval. 

This conhrms, hrst, that our integration procedure correctly hnds the roots of the equation 
Aoii = 0, and, second, that the critical exponent a = 7/6 can be achieved with the htting 
boundary conditions only, while in our case of condensate formation f'{0) = 0. This can 
explain also why the hnal distribution observed in Ref. was not resolved from a = 7/6. 
Indeed, the minimum value of the energy on the grid in Ref. was only £min = 10“^, while 
we have £min = 10“®, so the results of Ref. are more effected by the boundary conditions 


at small e and this happens in a shorter time. 


3. From the analytical point of view the 7/6 critical exponent corresponds to a stationary 
distribution, i.e. to a solution of the equation Jcoii = 0, while the shape of the self-similar 
distribution is dehned by Eq. ( |5.3b|) below, with yet to be determined parameter (3 which 
does not have to coincide with 7/6. 

The time dependence of the distribution function at e = 0 is shown in Fig. by the 
solid line. Vertical dot-dashed line corresponds to r = Tc and the distribution function 
asymptotically tend to this line according to the law /(O, r) oc (tc — t)“^'®, which is shown 
by the dotted line on this hgure. It is possible to hnd this time dependence at £ = 0 
analytically, using the self-similarity of the solution, see Section 


D. Condensation in the presence of condensate 

In the previous subsection we have seen that the distribution function of particles prior to 
the condensate formation tend to the power law with the exponent a been almost constant, 
see Fig. As one set of initial conditions for the condensation in the presence of the 
condensate we took the hnal distribution shown in Fig. |^. Though the limiting value of 
7/6 was not reached prior to condensate formation, this exponent is the special one and 
represents the root of the equation Ip = 0, see Appendix p. Because of that we have chosen 
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the function 


(3.6) 


/(e) =Cir-"''e-' 


for the second set of the initial conditions while integrating Eqs.( |2.9|) . In Eq. (|3.6|) C is 
some normalization constant, C = 1 corresponds roughly to the magnitude of the final 
distribution at £ ~ £min presented in Fig. ||a. We took initially Uc -C Utot- We will see 
shortly that the particular choice for ric, as far as this condition is satished, is insignificant. 


First let ns discuss the case with the initial condition (|3.6|) . The results of the numerical 
integration of the system (|2.9|) are presented in Fig. ^ The dashed line corresponds to 
the initial distribution. Solid lines correspond to the distribution function f{e) at different 
moments of time. Basically, evolution proceeds in the following way. First, the power 
law f{e) oc changes to the law f{e) oc 1 /e at small energies. And then this change 

propagates to the region of larger energies, see Figs. 9,10. Later on the power law stays 
at the equilibrium value a = 1, but the amplitude of the distribution function gradually 
decreases. The same kind of behavior of f{e,t) was found also in the paper 
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Again, we see that the essential part of the curves in Fig. ||b repeats itself under transla¬ 
tion from left to right and the evolution is self-similar. During this epoch (before “the wave 
of change” had reached the exponential tail of the initial distribution at £ > 1) approxi¬ 
mately 40% of particles had condensed. And, what is important, the number of particles in 
the condensate grows linearly with time at this epoch, ric/ntot = Bt. We found B ^ 2C^. 
This enables ns to eliminate the ambiguity in the initial value for Uc since B does not depend 
upon it. Indeed, in our simulations which were done in a hnite energy interval, during the 
hrst several iterations the system self-adjusts: A proper prohle of the distribution forms 
while the condensate reaches a particular value of ric- We can disregard this period and 
extrapolate the curves in Fig. ||b, ndr), and the self-similar character of the evolution back 
in time and to the region of smaller energies. 


The fraction of particles in condensate as a function of time is presented in Fig. ^ 

With the hnal distribution in Fig. ^ taken as the initial condition for the condensation 
in the presence of the condensate the fraction of condensed particles increases with time 
dramatically differently initially. Instead of ric/ntot oc r we obtain now ric/ntot oc ^/T. Both 
exponents can be understood analyzing self-similar solutions, see Section [VB| . 


IV. INTERACTION WITH ZERO TEMPERATURE FERMION BATH 


In this Section we consider ideal bose gas coupled to a fermion heat bath of inhnite 
capacity. The function F in the kinetic equation (|2.6|) in this case is dehned by Eq. ( 2.2b ) 
and 


For simplicity and dehniteness we shall consider fermions at zero temperature, Tbath = 0. 
The corresponding distribution function of fermions has the form ri{e) = — e), where /i 

is the Fermi energy. Such a simple form of the fermion distribution function allows ns to 
make one more analytical integration in the collision integral; at the end only one integra¬ 
tion remains to be calculated numerically. The price being paid for that is the increased 
complexity of the final analytical expression which we put therefore in Appendix p. 
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As an initial distribntion we choose the fnnction which has an eqnilibrinm form: 


finis) = 


exp [(£ - fiin)/Tin] - 1 ’ 


(4.1) 


where Tin and pin are initial temperatnre and chemical potential of the Bose system corre¬ 
spondingly. In what follows we choose the energy scale Eq and dehne the particle energy, the 
chemical potential and the temperatnre in nnits of Eq. To this end we choose /in(£ = 1) = 1) 
which gives ns the relation pin = [1 —log(2)]Tin. Under this condition only one dimensionless 
parameter Tin is left to dehne the initial state of Bose particles. 


The parameter Eq after rescaling will not enter the kinetic eqnation explicitly, bnt will 
dehne the time scale of the evolntion. We dehne the dimensionless time r as 


T=^^t 

64 : 71^771 


(4.2) 


We had specihed the initial distribntion Eq. ( |4.1| ) choosing Tin = 2 . For the Fermi 
energy we have nsed fip = 20. Since Tbath = 0, all Bose particles has to condense eventnally. 


We integrated kinetic eqnation in the energy interval 10“® < e < 100. We dehned the 
distribntion fnnction on the grid of 661 points eqnally spaced in the logarithm of energy and 
nsed Spline interpolation when calcnlating distribution function at intermediate points. We 
had implied Simpson algorithm on the grid of the same size for the collision integral . The 
non-conservation of particle number over the integration time t ^ 20 was Sn/n = 0.001%. 


The results of numerical integration of the kinetic equation (|C1|) are presented on Fig. 
r0| , where we plot the distribution function at several diherent moments of time. The 
distribution function behaves completely differently now as compared to the case of self¬ 
interacting bosons. It narrows in time gradually approaching a 6 - function 


The value of the distribution function at zero energy (see Fig. tend to inhnity with 
time only as r^. Consequently, it would require inhnite time for the distribution function 
to reach the inhnity. This does not necessarily means that the condensation time is inhnite 
in this case (see the corresponding discussion in the Introduction), but means that in this 
simple approach we can not hnd the condensation time-scale for the ideal Bose gas coupled 
to the thermal bath. However, we still can hnd the upper limit on the condensation time- 
scale, even without the knowledge when the phase transition to the condensed phase takes 
place actually. Indeed, the interaction with fermions necessarily leads to the self-interaction 
in higher orders of the perturbation theory (consider the box diagram constructed from 
the diagram Fig. ^d) with A ~ When the product g^fiO) exceeds 0(1), this induced 
self-interaction starts to dominate, and the problem reduces to the previous case of self¬ 
interacting bosons with the explosive grows of the distribution function. 


V. SELF-SIMILARITY AND ANALYTICAL SOLUTIONS 

The results of the numerical integration presented in previous sections suggest that during 
condensate formation the distribution function evolve in a self-similar way. Assuming self¬ 
similarity from the very beginning we can simplify the kinetic equation and then a number 
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of analytical results can be obtained In this section we discuss this analytical 


approach and compare it to the numerical results. 
We parametrize self-similar solutions as 


/(e.r) = A ‘^{r)f,(l) , 


( 6 . 1 ) 


where e = e/A{ t) and (3 is some constant. For dehniteness we shall assume /? > 0, i.e. in the 
following subsections A and C, where we shall consider the epoch prior to the condensate 
formation (when /(0,r) has to grow), A{t) —>■ 0 as time r increases , while A{t) will be 
growing function of time in subsection B which describes the epoch after the condensate 
has formed. We always can choose the normalization /s(0) = 1 if /^(O) is hnite, which is the 
case in subsections A and C, and A(0) = 1 otherwise. 

Self-similar solutions appear when the phase-space density is large, f + 1 ^ f, m which 
case the kinetic equation became invariant with respect to the rescaling. This condition is 
satished around the condensate formation, and as we have seen, the system approaches self¬ 
similar solutions starting with arbitrary initial conditions. We assume /(e) 3> 1 everywhere 
below. 


A. Self-interacting bosons prior to condensate formation 


Substituting parametrization (^) into kinetic equation (p.6|) we obtain: 




Separating variables we hnd 


_ 4 2 / 3-3 4 _ ^ 

W,(E) H- 


de 


C. 


(5.2) 


(5.3a) 

(5.3b) 


where Cg is (positive) separation constant. Integrating the hrst equation in (|5.3|) we obtain: 


3(T) = |2C.(T,-r)(/3-l)]‘''2('’-') , 


(5.4) 


where Tc is the integration constant apparently corresponding to the moment of time of the 
condensate formation. Indeed, the distribution function at zero momentum as a function of 
time is /(O, r) = A(r)“^, and /(O, r) has the pole at r = Tc if /3 > 1. 

We can conclude from Fig. ^a that at the tail where the distribution function is a power 
law, / cx it does not depend upon time. According to Eq. (0) this means that jS = a. 
This concludes the derivation of the function /(O, r) for the self-similar solution, which is 
plotted in Figs. 0(a,b) by the dotted line using a = 1.24. This value of a can be extracted 
from Fig. |b . The agreement of this function with direct numerical result is very good, 
see Fig. 0b. We conclude, that after the solution have reached self-similar form, the time 
dependence of the distribution function at zero momentum is given by 
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/(O, r) oc (r, - oc (r, - t)-"* ® 


(5.5) 


which reaches inhnity at a hnite time. Note that with the critical value a = 7/6 we would 
obtain /(0,r) oc (tc — This function is plotted by the dashed line in Fig. and is 

quite different from what we observe numerically. 


B. Self-interacting bosons and the condensate 


Substituting parametrization Eq. (^TT|) into Eqs. (^.9| ) we obtain: 

= „^(r)A(Tr^»h . 

CLT 


(5.6a) 


A{t)A{t\ 




/3/s(di) + El 


dfs{£i) 

dii 


= -A{Tf ^^Ip{ei) - nc(r)v4(r)(^ , (5.6b) 


where Ic and Ipc are integrals in r.g.s. of Eqs. (p.9a)) and (|2.9b|) respectively. Under the 
condition 


n,(r) = 


variables in Eq. ( |5.6b| ) separates: 


i(r)7l(r)-^-' = A{Tf-^l^Cs , 


(5.7) 


(5.8a) 


-a 


Pfs{.£i) + £i 


dfsi£i) 
ddi 


dpi^l) T C'^/pc(£'i) 


(5.8b) 


If we solve Eq. ( b.TD with respect to A{t) and substitute it into Eq. (|5.6a|) we hnd 

(3-2/3)/4(/3-l) 


nrir) = C„ 


4(l-/3)/c 
(2/3 - 3) 


■ r 


(5.9) 


which gives the time dependence of the number of condensed particles. As a self-consistency 
check we hnd that the same substitution, but with Eq. (|5.7|) solved with respect to ndj), 


gives the condition equivalent to Eq. (|5.8a]) if Cg = 2Ic/{3 — 2/3). This dehnes Cg- The 


solution of Eq. (|5.8a|) is: 


A(t) = 


. , 4(1 - P)h ' 
( 2 / 3 - 3 ) 


l/2(/3-l) 


(5.10) 


As we have seen, the equilibrium prohle f oc. 1/e propagates from the region of small 
energies to the region of large energies as a “wave”, leaving / cx to be time independent 
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in the region which this “wave” has not reached yet. This, again, means that (5 = a. With 
a = 7/6 as an initial condition we hnd from Eq. ( ^.91) that the condensate linearly grows with 
time, nc(r) = const x r. This is in agreement with the results of our numerical experiment, 
see Fig. ^ This self-similar regime persisted till the “wave” have reached exponential tail 
where a is different. During this period ~ 40% of particles had condensed. With the initial 
condition a = 1.24 we obtain dramatically different rate of condensate growth, nc(r) oc 
which, again, was conhrmed by the direct numerical integration. 


C. Interaction with degenerate Fermi gas 


Now we substitute parametrization (ED into the kinetic equation with the function 
F given by Eq. ( 2.2b ). Left-hand side of equation (^) remains unchanged in this case, 
however, its right-hand side reads now This changes equation (|5.3| a) to 


A = -A^-^ a 


(5.11) 


In the case of the self-interacting bosons we were able to hnd (3 using the information gained 
from the numerical integration that on the tail of the self-similar distribution the value of / 
does not depend upon time. That gave us [3 = a and we have extracted a as a logarithmic 
derivative of the distribution function with respect to the energy. We could not use integrals 
of the motion in that case since the particle number (as well as the energy density) saturated 
at the very end of the distribution where the self-similarity does not hold anymore (since 
on the self-similar solution we had a < 3/2). Now, on the contrary, the particle number 
saturates already near the core of the self-similar distribution, see Fig. (moreover, a 
grows with time). This enables us to hnd (3 using conservation of particles. 

Substituting parametrization (|5.1|) into Eq. (|3.4b|) we hnd (3 = 3/2. Consequently, the 
equation (|5.11|) reads now A = C^. We hnd hnally that 


f(e = 0,T) = A oc (t + T„)= 


(5.12) 


This is in good agreement with the results of our numerical integration, see Fig. |^. 

The ^-dependence of the distribution function can be found from the equation ( |5.3bD , 
which in the limit of large /i takes the form 


3/s 

2 ^ di 


A 


1 

/ dxx^'‘^{x — e)fs{x) + / dx{x — e)fs{x) 

Ve Jo Je 


( 5 . 13 ) 


Let us make here the connection to results of Ref. . One can see that it is possible to 
rewrite the hnal expressions of [|^] in a self-similar form, Eq. (EID. with the same value of 
/? = 3/2 and the same function A{t) cx at large r. This gives fs = fs{\/sT). Using this 
and noting that Levich and Yakhot started from the initial conditions corresponding to the 
critical point, /s(d) cx 1/d, we hnd their result f{e, r) cx r. Our direct numerical integrations 
with an appropriate initial conditions conhrm that. However, in the runs the results of which 
we had reported above we have had diherent boundary conditions, /^(O) = 1, which leads 
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to /(O, r) oc at large r, Eq. ( ^■12| ). We see that the system which is father away from the 
critical point condenses with a higher rate. 


The solntion fonnd in Ref. for fg, which with onr initial and bonndary conditions 


corresponds to /s(d) = exp(—cr\/l), is not an exact solntion to eqnation ( p.RlI ), however, it 
is a reasonable approximation. This solntion is plotted in Fig. 0 by the dotted line. 


VI. CONCLUSIONS AND DISCUSSION 


We have stndied numerically the kinetics of condensation of the weakly interacting Bose 
gas. We have found that the distribution function evolves very differently in cases of Bose- 
bath interaction and self-interaction of Bose particles. This results from the different struc¬ 
ture of the kinetic equations describing particle-bath and particle-particle interactions [the 
latter contains extra powers of f{e)]. In the hrst case the distribution function of the excess 
particles, which eventually form the condensate, narrows gradually with time and continu¬ 
ously approaches a 6 - function in inhnite time, as it was found in Ref. assuming small 
energy exchange per collision. 


For the case of self-interacting bosons we hnd, in qualitative agreement with Refs. [25,16 


that the singular power law prohle f{e) oc of the distribution function forms in a 

hnite time. This is close to f{e) oc which is a stationary solution of the kinetic 

equation and corresponds to the constant flux of particles in momentum space towards the 
condensate. Different exponents result from different boundary conditions during the Bose 
condensation, namely df/de{0,t) = 0, as opposed to the case of stationary turbulence. 
While the exponents differ by only 6% this leads to signihcantly different values for other 
critical exponents, e.g. f3 in /(0,t) oc {tc — t)~^ , or for the rate of the condensate growth 
7 in Uc oc P. We have found f3 = 2.6 and 7 = 0.54, while for f{e) oc we would have 
(3 = 3.5 and 7 = 1 . 


One could expect that this very natural regime of the constant flow of particles into the 
condensate will persist in the presence of the condensate as well until all excess particles 
from the high energy tail are in the ground state. Nevertheless, this is not the case and 
at the moment the condensate appears, its presence terminates this regime. Instead of 
such a steady flow through the entire energy interval, particles from all energy levels jump 
directly to the condensate, while the remaining particles maintain an equilibrium shape of 
the distribution function f{e) oc The constant of proportionality in this law gradually 
decreases until it reaches an equilibrium value. 


The build up of coherence cannot be observed in the framework of the Boltzmann kinetic 
equation, but the kinetic description has to be valid prior to and after the moment of 
condensate formation. We have seen that evolution at both stages is self-similar. This 
allows US to obtain a number of useful analytical relations, e.g. the time dependence of the 
distribution function near the point of condensate formation, Eq. (|5.4|) , and Eq. (|5.9| ) for 
the subsequent growth of the condensate. We have found the duration of both stages, which 
is hnite and of order ~ 20. 


Our results of numerical integration of the Boltzmann equation might be used as an 
initial condition in a theory of condensate nucleation. 
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Another interesting application of our results could be in a description of Bose-star 
formation. A complete theory of this process requires the inclusion of gravity and the 
problem is more involved. While the isotropy in momentum space will be broken, the kinetics 
of Bose relaxation in the collapsing core can proceed largely along the lines described in the 
present paper. Moreover, the relaxation time can become shorter, since the spatial collapse 
of the core under the influence of self-gravity would cause an additional increase of the 
density of particles, which would speed up the relaxation. 


It is remarkable that in spite of the apparent smallness of the axion quartic self-coupling, 
A". ~ 10 (where /12 = /a/10^^ GeV) the relaxation time in axion miniclusters can 

be comparable to the age of the universe 0,|1. This is because the relevant quantity is not 
the self-coupling itself, but the product of A with the mean held strength squared (or with 
the mean phase-space density of “particles”, which in the present case would rather be the 
phase-space density of classical scalar waves). For particles bound in a gravitational well, it 
is convenient to rewrite the expression for the relaxation time in the form Tc ~ m7 
0, where p is the energy density in a minicluster and Ve is the escape velocity. For axions, 
the relaxation time is smaller than the present age of the Universe if the energy density 
satishes pio > lO^vs^/Ju, where pio = p/(10eV)^ and vs = Ue/10“®. It was shown in 
Refs. 0 that those large densities can be achieved naturally in axion miniclusters and, 
therefore, a “Bose-star” can form . 

Our results are applicable whenever Eq. ( p.l|) holds and might also be interesting for a 
study of, e.g., Kolmogorov turbulence. 
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APPENDIX A: REDUCTION OF THE COLLISION INTEGRAL 


In this Appendix we will make as many integrations in the collision integral analytically, 
as it is possible. In the case we are interested in, he., four-particle interaction vertex and 
the isotropic one-particle distribution function, seven out of nine integrations can be done 
and only two integrals upon the energies of incoming particles remain the for numerical 
treatment. We start from the non-relativistic counterpart of Eq. (p.l|) 


where F{f) is dehned in ( p.2|) with the distribution function / which depends upon particle 
energy £ only. We shall make use of the identity 
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d^X 


and explicitly separate out angle integrations 

d^Pi = d(f)id cos Oip'fdpi = mpidVLidSi . 
Collision integral (|A1|) takes the form: 


draW - 


\Mf^ 

Gdvr^m 


h(£l + £2 — — £2)F{f)Dd£[d£2d£2 ■ 


(A2) 


(A3) 


(A4) 


where we have assumed that the matrix element does not depend upon momenta and we 
have dehned 

D = J X^dX J J J . (A5) 


All angle integrals are trivial in (|A5|) and we get 


D = 


^ sin(Api) sin(Ap 2 ) sin(Api) sin(Ap 2 ) • 


(A6) 


vrpi Jo A^ 

After some algebra and taking into account the energy conservation condition this gives 

D = ^min[pi,p2,Pi,P2] • (A7) 

Now, the energy d-function can be trivially integrated away and we are left with two- 

(A8) 


dimensional integral upon energies of incoming particles 

d/(£i) 


dt 




Fif)Dd£[d£'2 , 


where £2 = £'i — £'2 — £i- 


APPENDIX B: STATIONARY SOLUTIONS TO THE KINETIC EQUATION 


1. Kolmogorov’s spectra 

In this Appendix we shall describe power law solutions of the equation Jcoii = 0. We shall 
consider the region of large values of the distribution function, /(e) 3> 1, he., F{f) will be 
given by Eq. (|3.5|) . The kinetic equation with this form of F{f) was extensively studied in 
the literature devoted to the problem of the Kolmogorov turbulence and in Ref. all 
power law solutions to the equation Jcoii = 0 were found. Two of those solutions are evident. 
One is just the limit of the Eq.( p.5|) at £ —> 0, and Uc = 0 and has the form / = const/e. 
Another one is / = const. Both are solutions to a simpler equation F{f)6{J2^) = 0. To 
hnd the remaining solutions is a non-trivial task and with the distribution function of the 
form /(e) = the trick is to map all labelled regions in Fig. ^ to the region I by means 
of Zakharov transformations [11^1: 
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II I 


III I 


IV ^ I 


^2 


£i4/£2 


^2 


^2^i/4 

- eVe', 


4 ^ ^?/4 
4 ^ ^2^iM 


where 82 = + £2 ~ ^i- Under this transformations the snm of integrals over regions I - IV 

can be represented as an integral over the region I only: 



(Bl) 


where 82 = + £2 ~ cind /3 = —7/2 + 3a. We see that in fact there are fonr solntions to 

the eqnation Icon = 0 which correspond to a = 0, 1 or /? = 0, 1. This gives 

a = 0, a = 1, a = 7/6, a = 3/2 . (B2) 

Onr nnmerical self-similar solntion, Fig. |^, tend to the power law / = £“^6. however, it 
never reaches this critical exponent prior to the condensate formation, see Figs. 7 and 8. 
Physically the solntion with a = 7/6 corresponds to the constant flnx of the particle nnmber 
from the region of lager energies towards the fnture condensate. 


2. Solutions in the presence of condensate 

In the presence of the condensate the nnmber and the natnre of zeros of the collision 
integral changes. The kinetic eqnation can be represented as a system of two eqnations, 
Eqs. (|2.9|) . The collision integral in Eq. (p.9b| ) consists of two terms. The hrst term. Ip, is 
the same as in the absence of the condensate, and the corresponding integral is zero on 
the solntions corresponding to the same valnes of a, Eq. (P^) . However, the second term 
does not. Assnming the power law /(e) = the integration in this term can be made 
analytically with the resnlt proportional to 

,i- 2 ar 4 4"F(l-a)(7r + F(|-a)F(a-|)) l 

[a-1 ^ 2v/?F(|-a) J ‘ ^ ^ 

This expression as a fnnction of a has two zeros only at 

a = 1 and a = 3/2 , (B4) 

and conseqnently the total collision integral in Eq.( p.9b| ) is zero on those exponents only. 
Moreover, the collision integral in Eq. (|2.9a|) is zero with a = 1 only. This explains why the 
power law /(e) = e“^/® changes to /(e) = e“^ in the presence of the condensate, see Fig. 

lb. 
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APPENDIX C: KINETIC EQUATION FOR IDEAL BOSONS INTERACTING 

WITH DEGENERATE FERMI GAS 


Assuming fermions to be at zero temperature, be., their distribution function to be given 
by v{^) = ~ allows us to make one more integration in Eq. (|2.6|) analytically. After 

some algebra we obtain: 

+ J 2 ) + [/(^) + 1]('^3 + J 4 .), (Cl) 


where 

Ji = — f dxx^^'^{e — x)[f{x) + 1] , (C2a) 

we Jo 


^ {Jo + e- fi) + l]- j dx[f{x) + 1] 

- J^ dxx^/^{fi-^)[f{x) + l]j (C2b) 

J 3 = r]{e) {^J^ dx{x - e)f{x) + J dx{iJ - ^)f{x) - dxx^^^ f{iJ, + e - x)^ (C 2 c) 

J 4 = —^ Iy dx/j,^^"^f{x) — J dxx^''"^f{iJ, + e — x)'^ (C 2 d) 

where /i is the Fermi energy. This equation was subject to the numerical integration in 
Section 
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FIG. 1. The scattering processes 

we consider 

in Eq. (2.1). (a) is the case of self-interacting 


bosons, and (b) is the case of ideal Bose gas interacting with the cold fermion bath. 



FIG. 2. The white areas correspond to the integration region in the collision integral, Eq. (2.6) 



FIG. 3. The phase plane of the system. Below the solid line /x = 0 the condensate is non-zero 
in equilibrium. The dashed line corresponds to the variation of the parameter /q. 
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FIG. 4. Snapshots of the distribution function in non-condensate case. The dashed line cor¬ 
responds to the initial distribution; the equilibrium distribution is presented by the solid curve; 
dotted curves correspond to the results of numerical integration at different moments of time 
r = 1,5,10,20. 



FIG. 5. The distribution function at zero energy as a function of time for two cases presented 
in Fig. 4. 
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FIG. 6. The distribution function of self-interacting bosons (top) and its logarithmic derivative 
(bottom) are shown at different moments of time prior to condensate formation. The dashed line 
is the initial distribution; the dotted line corresponds to the power law /(e) oc 


26 




8 10 12 14 16 18 20 



0.01 


0.1 1 

('Tc -t) 


FIG. 7. The distribution function at e = 0 for the case presented in Fig. 6 (the solid curves) 
The dotted curve in both figures is a solution to self-similarity equations, / oc (tc — whicl 

was found with a = 1.24, see Section V. The dashed line in Fig. (b) is the power law / oc 
which would correspond to a = 7/6. 
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FIG. 9. Fraction of condensed particles as a function of time for the case presented in Fig. 8. 



FIG. 10. The distribution function of bosons interacting with the cold gas of fermions is shown 
at different moments of time by the solid curves. The dashed curve is the initial distribution, the 
dotted curve corresponds to the Fokker-Planck approximation of Ref. [13]. 
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